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ABSTRACT 

We report two-dimensional MHD simulations which demonstrate that pho- 
tospheric 5-min oscillations can leak into the chromosphere inside small-scale 
vertical magnetic flux tubes. The results of our numerical experiments are com- 
patible with those inferred from simultaneous spectropolarimetric observations 
of the photosphere and chromosphere obtained with the Tenerife Infrared Po- 
larimeter (TIP) at 10830 A. We conclude that the efficiency of energy exchange 
by radiation in the solar photosphere can lead to a significant reduction of the 
cut-off frequency and may allow for the propagation of the 5 minutes waves ver- 
tically into the chromosphere. 

Subject headings: MHD; Sun: magnetic fields; Sun: oscillations; Sun: photo- 
sphere, Sun: chromosphere 



1. Introduction 

Several recent investigations ha ye reported the detection of 5 minute oscillations in the 
chromosphere and corona of the Sun (Kriiger et aliboOl : De Moortel et al. 2002 : De Pontieu et al 
20031 : ICenteno et al~1l2006al : IVecchio et alJl2007h . These 5-min oscillations have been detected 
mainly above facular and network areas, whil e oscillations in the outer atmospheric regions 



of sunspots show mainly a 3-min periodicity (ICenteno et al.ll2006bl ; iBogdan fc Judge! 12006 
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and references therein) , similar to what it is found in the internetwork regions of the Sun. An 
ex planation for the shift in the wave period with height from 5 to 3 minutes was suggested 
by iFleck fc Schmitzl (Il99ll ). who argue that this is a basic phenomenon due to resonant 
excitation at the atmospheric cut-off frequency (see also Fleck & Schmitz 1993). The low 
temperatures of the high photosphere give rise to the 3-minute cut-off period. However, it is 
unclear why such a shift in the wave period is not observed in small-scale structures present 
in facular and network regions. How can the evanescent 5 minute oscillations propagate up 
to the chromospheric heights? 



De Pontieu et al.l (120041 ) argue that the inclination of the magnetic field is essential for 



the leakage of p-modes strong enough t o produce the dynam ic jets observed in active re- 
gion fibrils. Following the same concept. I Jefferies et al.l (120061 ) claim that inclined flux tubes 
explain their lower chromospheric observations of propagating waves (whatever their ampli- 
tude). If (acoustic or slow MHD) waves have a preferred direction of propagation defined by 
the magnetic field, the effective cut-off frequency is lowered by the cosine of the inclination 
angle with respect to the solar local vertical. This allows evanescent waves to propagate. 
However, there are also works reporting vertically propagating 5 minute waves in facular 
and network regions at chromo spheric heights that cannot be easily exp l ained by the above- 
mentioned mechanism (see e.g-. lKrijger et al.ll200ll ; ICenteno et al.ll2006al ; IVecchio et al.ll2007l . 
and references therein). Alternatively, a decrease in the effective acoustic cut-off frequency 
can be produced by taking into acco unt the radiat i ve energy losses, with out the need of 
assuming inclined magnetic flux tubes (|Robertslll983l ; ICenteno et al.ll2006al ). If the radiative 
relaxation time (trr) is sufficiently small (as expected for small-scale magnetic structures), 
the cut-off frequency can be reduced allowing the 5-min evanescent waves to propagate. 



In this letter, we consider the obser vational eviden ces provided by lCenteno et al.l (j2006al ). 
and extend the theoretical analysis by iRobertsi (119831 ) with the help of non-adiabatic, non- 
linear 2D numerical simulations of magneto-acoustic waves in small-scale flux tubes with a 
realistic magnetic field configuration. 



2. Summary of the spectropolarimetric observations 



A facular region, located close to the disc center (// = 0.95), was observed at the German 
Vacuum Tower Telescope (VT T) of the Observatorio del Teide on 14 July 2002, using the 
Tenerife Infrared Polarimeter (IMartmez Pillet et al.lll999l ). The data set consists of a time 
series of 81 minutes (with a temporal cadence of 5.4 s) taken with a fixed slit position. The 
slit was 0.5 wide and 40" long. The full Stokes profiles were measured in a spectral range 
that spanned from 10825.5 to 10833 A, with a spectral sampling of 31 mA per pixel. This 
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spectral region includes a photospheric Si I line at 10827.09 A and a chromospheric Helium 
I 10830 A line, which is indeed a triplet whose blue component (A 10829. 09A) is quite 
weak and difficult to see in an intensity spectrum, and whose red components (A 10830.25, 
A 10830.34 A) appear blended. Several non-LTE radiative transfer investigations indicate 
that the He I 10830 A spectral line radiation is originated in a relatively thin layer in the 



upper chromosphe r e, about 2000 km abo ve the base of the photosphere (jPozhaloval 11988 



Avrett et al.l Il994l ; ICenteno et al.l l2008bl ) . The LOS velocity variations of these spectral 
lines have been the main object of the analysis. They were retrieved by applying a Milne- 
Eddington (ME) inversion to the chromospheric Helium lines and a standard LTE inversion 
code to the photospheric Sil i con line. For a detailed description of the observations and their 
analysis, see ICenteno et al.l (j2006al . l2008al ). Below we summarize briefly the main findings. 



The left panel of Fig. [T] shows the power spectra of velocity averaged over all the posi- 
tions along the slit in the photosphere (dashed line) and the chromosphere (solid line) as a 
function of the frequency of the oscillating mode. The right panel represents the difference 
in phase (S(j)) between the chromospheric and the photospheric velocity oscillations, also as 
a function of frequency. It is interesting to note that the photospheric and chromospheric 
facular power spectra show very similar features with peaks at the same positions, having 
their main contribution in the 3 mHz band (5- minute oscillations). The facular phase differ- 
ence spectrum shows a very noisy behavior below 2 mHz, indicating that there is no wave 
propagation in this frequency regime. Above this point (and also around the 3 mHz band) 
the phase difference starts to increase with frequency, meaning that these frequency modes 
do propagate from the photosphere, reaching the chromosphere some time later. 



Theoretical considerations 



The low amplitudes of the observed linear polarization signals of the Si I line (not larger 
than 2 x 10~ 3 in units of the continuum intensity) imply inclination angles smaller than 
10-15 degrees to the solar vertical. Consequently, throughout the present paper we make the 
assumption that the magnetic field observed in the facular region is close to the local vertical. 
The equations describing the vertical propagation of a wave in a vertical magnetic field in 
a gravitationally stratified isothe rmal atmosphere are formally equal to th ose governing the 
propagation of an acoustic wave (IFerraro fc Plumptonlll958l ; lRobertsll2006l ). These equations 
have the following well known solution describing the variation of the wave amplitude V(z) 
with height: 

V(z) = V e z/2H e k > z , (1) 
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Fig. I. — Left: Photospheric (solid line) and chromospheric (dashed line) LOS velocity 
power spectra in the observed plage region. Right: Phase spectrum for the LOS velocity 
oscillations. Each cross on the figure is computed as the difference from the chromospheric 
phase and the photospheric phase for one frequency and one position along the slit. Solid 
lines represent the best fit of the theoretical model with values T=9500 K, t rr = 10 s, and 
Sz = 1500 km. 
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Fig. 2. — Phase difference between waves at different heights calculated theoretically as 
a function of frequency. Dashed line: assuming adiabatic oscillations; solid line: using a 
radiative relaxation time of trr = 10 s. Note the difference in the effective cut-off frequency 
between the two cases. The parameters of the calculation are 5z = 1500 km, T =9000 K. 
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where kl 



- u?)/c 2 cl = i*gH, ul = i*g/4H and 7* = (1 - ijUT RR )/(l - iur RR ) 



( IMihalas &: Mihalaall984j ). In general, the vertical wave vector k z is a complex number. If 
the wave propagation is adiabatic (t rr — > 00), the real part of the wave vector k r becomes 
equal to zero for frequencies below the cut-off and no propagation is possible. However, if 
energy exchange by radiat ion is taken in to account -that is, if the radiative relaxation time 
t rr has a finite value (see ISpiegellll957l ). k z always has both an imaginary and a real part. 
Thus, formally, the wave propagation is possible at all frequencies in this case. It should 
be noted that the rad iative relaxation time is very small in small-scale structures of the 
solar photosphere (see iKneer fc Trujillo-Buenol 119871 ). The properties of the wave solution 
are illustrated in Fig. [2j It shows the phase difference 5<p = k r 5z between the waves at 
two heights (separated 1500 km) as a function of frequency for the cases r RR — > 00 (dashed 
line) and t rr —10 s (solid line). An effective cut-off frequency may be defined for the case of 
t rr = 10 s, being significantly lower than for the adiabatic case. This fact causes important 
effects on wave propagation. 

The solid curve on the right panel of Fig. [1] gives the best fit of the above-mentioned 
model to the observed phase spectrum in the facular region. The free parameters of the 
model are the temperature, the radiative relaxation time and the height difference between 
the two velocity signals. The best fit results in a value of t rr = 10 s. 



MHD simulations 



To generalize these conclusions (based on a simple model of linear wave propagation) 
to the case of a more realistic atmosphere with horizontal and vertical stratification in 
all the physical parameters, a complex magnetic field configuration and non-linear wave 
propagation, we solved numerically the set of MHD equations in this realistic environ- 
ment. The numerical code used fo r the solution of these equations is described elsewhere 
( IKhomenko fc Colladoa 120061 120071 ). Radiative losses were taken into account by means of 
Newton's law of cooling. The initial magn e tostat ic flux tube model was constructed following 
the method described by iPneuman et al.l (119861 ) (see Fig. [3]). 



Previous investigations showed that the photospheric driver that excites oscillations 
inside the flux tube should have a vertical component in orde r to generate an acoustic fast 
wave already in the photosphere (IKhomenko fc Colladosl 120071 ). If this were not the case, no 
correlation would be observed between the photospheric and chromospheric velocity signals, 
contrarily to what is suggested by observations (see the power spectra of Fig. [[]). 



Thus, we applied a vertical periodic driver at the bottom boundary of the simulation 
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Fig. 3. — Magnetostatic flux tube model. The background image quantifies the temperature 
values. The magnetic field along the flux tube axis varies with height from 740 G in the 
photosphere to 37 G at z=2000 km, the plasma f3 decreases from 4 to 7xl0 -4 and the cut-off 
period increases from 230 s in the photosphere to 280 s in the chromosphere. 
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domain by varying the vertical velocity as: V z (x, t) = Vosin(ut) x exp( — (x — xo) 2 / 2a 2 ), 
where a = 160 km is the horizontal size of the pulse, Vq = 200 m s _1 is the initial amplitude 
at the photosphere and x corresponds to the position of the tube axis. The period of the 
driver was 300 s (i.e., it is above the cut off period). In what follows we will compare 
two simulation runs that are identical except for the value of the radiative relaxation time. 
While the first run is in the adiabatic regime (trr — > oo), the second one was carried out 
with trr = 10 s (constant through the whole atmosphere). 

Fig. H] gives time-distance plots of the vertical velocity obtained from the simulations 
at two heights: the upper photosphere around 400 km (top) and the chromosphere around 
1500 km (bottom). The vertical photospheric driver generates a fast magneto-acoustic wave 
that is essentially acoustic in the photosphere where (3 > 1. This wave propagates upwards 
through the (3 = 1 transition layer, preserving its acoustic nature and being transformed into 
a slow magneto-acoustic mode higher up. 

In the adiabatic case (left panels of Fig. 4), we can see that the velocity in the photo- 
sphere has a well defined 5-minute periodicity within the flux tube (horizontal distance from 
-200 to 200 km). Higher up in the chromosphere, this behavior changes. Note that it takes 
some time for the oscillations to reach the chromospheric heights (bottom panels of Fig. H|). 
In the chromosphere the tube expands and occupies now a larger area in the time-distance 
plot (from about -400 to 400 km). At the beginning of the time series, two strong shock 
waves develop, and are followed in time by strongly damped oscillations. Only weak shocks 
exist in the stationary phase of the simulations and the wave periodicity is 3 minutes. Thus, 
in the adiabatic case, there is a shift with height of the dominant period from five minutes 
in the photosphere to three minutes in the chromosphere. 

This behavior contrasts with that of the simulations with trr = 10 s (right panels 
of Fig. H]). Contrarily to what happens in the adiabatic case, the 5-minute periodicity is 
preserved both in the photosphere and in the chromosphere. Waves are linear in the deeper 
layers, they propagate vertically upwards through the temperature minimum and form shocks 
with 5-6 km s _1 amplitudes in the higher layers. 

The left panel of Fig. [5] shows the power spectra of the photospheric and chromospheric 
vertical velocity oscillations in the simulation run with trr = 10 s. Now there is no shift in 
the wave period as we go to the higher atmospheric layers, and both the photospheric and 
chromospheric power spectra are rather similar (right panel). This figure can be directly 
compared to the observed power spectra in Fig. [T] (left panel). The similarity between the 
observed and the simulated power spectra is notable. The right panel of Fig. [5] gives the 
difference in phase between the chromospheric and the photospheric velocity oscillations as 
a function of frequency, as calculated from the simulations. This figure can be compared 
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Fig. 4. — Time-distance plots of the vertical velocity from the simulations at photospheric 
(top) and chromospheric (bottom) heights. Left: adiabatic calculation; right: calculation 
with a radiative relaxation time trr = 10 s. Note how the horizontal extent of the tube is 
larger in the chromosphere because of the field expansion. The vertical bars on the right of 
each plot represent the velocity scale in km s _1 . 
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Fig. 5. — Left: photospheric (dashed line) and chromospheric (solid line) power spectra 
obtained from the simulation with trr equal to 10 s. Right: phase difference spectrum for 
the vertical velocity oscillations. Each cross on the figure is computed as the difference from 
the chromospheric phase and the photospheric phase for one frequency and one position 
inside the flux tube. 
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with the observed phase spectra in Fig. [T] (right panel). It can be seen that the simulations 
describe correctly the phase speed of waves in the frequency range 0-5 mHz, as well as the 
value of the effective cut-off frequency. This gives additional support to our conclusions. 
Radiative losses are expected to play an important role in small-scale magnetic structures, 
such as those present in facular regions and, as already shown, they are able to decrease the 
cut-off frequency allowing the five-minute oscillation to propagate into the chromosphere. 

The results shown above were obtained with a fixed value of the radiative relaxation time 
{trr — 10 s). We have checked that similar conclusions are reached with a height-dependent 
relaxation time changing from 10 s in the photosphere to 100 s in the chromosphere. 

5. Conclusions 

In this paper we have performed an observational, analytical and numerical analysis 
of the propagation of the 5-min oscillations from the photosphere to the chromosphere in 
small-scale vertical magnetic flux tubes. From the observational point of view, it has been 
demonstrated that the power spectra of oscillations in the photosphere and the chromosphere 
in a facular region are similar and peak at 3 mHz. The observed velocity variations are 
compatible with the longitudinal acoustic wave propagation in a vertical magnetic field from 
the photosphere to the chromosphere. The analytical theory of acoustic waves in a stratified 
isothermal atmosphere embedded in a vertical magnetic field suggests that the effective cut- 
off frequency can be significantly reduced if the radiative losses of the oscillations are taken 
into account. This allows the propagation of the evanescent waves at 3 mHz. 

The numerical solution of the complete set of the MHD equations allows to generalize 
this conclusion to the situation of a complex magnetic field configuration, stratification in 
all atmospheric parameters and non-linear wave effects. We conclude that our photospheric 
and chromospheric observations can be explained assuming that the radiative relaxation 
time is sufficiently small (as expected in small-scale magnetic structures). This significantly 
reduces the cut-off frequency and makes evanescent 5-minute waves propagate vertically from 
the photosphere towards the chromosphere without changing their period with height. Our 
conclusion that vertical magnetic field concentrations can also channel the 5-min oscillatory 
power of the photosphere towards the chromosphere in facular regions extends the actual 
possibilities of the dynamic and magnetic couplings between these two important layers of 
the solar atmosphere. Additional work on three-dimensional MHD simulations with non- 
LTE radiative transfer physics will be necessary to investigate the role of this mechanism in 
more realistic models of the highly dynamic conditions of the magnetized solar atmosphere. 
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